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Abstract 

> 

. Inspired by Kitaev's argument that physical error correction is possible in a system 

of interacting anyons, we demonstrate that such "self-correction" is fairly common 
in spin systems with classical Hamiltonians that admit the Peierls argument and 
where errors are modelled by quantum perturbations. 
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1 Introduction 

Error correction is a key ingredient in any good recipe for quantum information 
processing (see [1] for an overview and further references). Many ingenious 
schemes have been invented to that effect. A particularly interesting approach 
has been suggested by Kitaev and colleagues in a series of beautiful papers 
[2-4], namely the possibility of implementing quantum error correction on the 
physical level. 

The idea is quite simple. In [2], Kitaev considers a specific Hamiltonian H\ for 
a quantum spin-1/2 system on a finite lattice AcZ 2 with periodic boundary 
conditions (the spins are taken to reside on the bonds of A). For any finite 
A, the ground state of H\ is fourfold degenerate; information is stored in the 
corresponding four- dimensional eigenspace. Addition of a small local pertur- 
bation of a certain kind modifies the Hamiltonian to H\(e), where e is the 
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perturbation strength. The effect of the perturbation is to introduce an en- 
ergy splitting between the degenerate ground-state levels of the unperturbed 
Hamiltonian. Kitaev then argues that there exists a constant e so that, at 
low temperatures and for all |e| < e , the energy splitting for sufficiently large 
A is given by exp (— cy/N^), where N\ is the number of bonds in A. In other 
words, in the thermodynamic limit A j Z 2 , the ground state is still fourfold 
degenerate, and any sufficiently weak perturbation is "washed out" by the 
system itself. 

The four-body interactions comprising Kitaev's Hamiltonian were originally 
considered by him in Ref. [5] as generators for a family of stabilizer codes [1], 
which he termed "toric codes." The remarkable feature of toric codes is the fact 
that, despite their apparent nonoptimality (in the sense of [6]), they require 
only local operations for their implementation and can correct any number of 
errors (by making the lattice large enough). The bulk of Kitaev's analysis of 
toric codes was concerned with their properties as "conventional" quantum 
error-correcting codes [7] that require active intervention through frequent 
measurements and other external processing. The issue of constructing "self- 
correcting" quantum spin systems on the basis of toric codes has been taken 
up again only very recently by Dennis, Kitaev, Landahl, and Preskill [4]. Their 
approach, however, is centered around topological features of toric codes and 
delves deep into such subjects as nonabelian gauge theory [4,8]. 

On the other hand, the very idea of physical error correction is so tantaliz- 
ing, both practically and conceptually, that one cannot help but wonder: how 
generic are phenomena of this kind? In this Letter I show that a few results in 
statistical mechanics of quantum spin systems point towards the conclusion 
that physical error correction is fairly common in such systems, under quite 
reasonable conditions. These preliminary findings are, in my opinion, sugges- 
tive and interesting enough, and might spur further research into this topic 
both in the quantum information community and in the statistical mechanics 
community. A detailed analysis is forthcoming [9]. 



2 Laying out the ingredients 

First of all, let us agree on the ingredients necessary for analysis of a self- 
correcting quantum spin system. Let A C Z" be a finite lattice. We assume 
here that v > 2 is an integer. Let H be the (2S + l)-dimensional Hilbert 
space of a single particle of spin S. Spins are situated on lattice sites £ G A 
(in Kitaev's construction, spins were located on lattice bonds). In order to 
retain a superficial analogy with stabilizer codes, we will assume that the 
unperturbed Hamiltonian is classical, i.e., the interactions comprising it 
generate an abelian subalgebra of the algebra B(Ha) of all linear operators on 
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Ha '■= ®£eA^j where He is an isomorphic copy of Ho- That is, 



(1) 



McA 



where each $ M is a Hermitian operator on 7i M := ® eeM Hi, and [<E»m, ^tv] = 
0. We assume periodic boundary conditions (that is, the lattice A is drawn 
on the torus Z"/Z). We let {\g_)} be the orthonormal basis of H\ in which 
H\ is diagonal; the basis vectors are labelled by classical spin configurations, 
a = {a(£)}eeA with a(£) G {—S, —S+l, . . . , S—l, S}. We also assume that the 
smallest eigenvalue of H\ is equal to zero, and that its geometric multiplicity 
is m > 2. We denote the corresponding eigenspace by H\. 

The effect of errors is modelled by introducing an off-diagonal perturbation 
term to the Hamiltonian: 



where e is a positive constant and P is a Hermitian operator whose exact form 
is, for the moment, left unspecified. Addition of the eP term will perturb the 
eigenvalues of H\, resulting in energy splitting between orthogonal ground 
states of the original (unperturbed) Hamiltonian. Consequently, we define 



The basic idea behind a self-correcting quantum spin system thus boils down 
to the following. Information is stored in the ground-state eigenspace of the 
unperturbed Hamiltonian H\. The multiplicity m is, obviously, dictated by the 
desired storage capacity: when m = 2 k , our "ground-state memory cell" will 
hold k qubits. Errors will cause some of the information to leak out into excited 
states. In order for error correction to take place, the system should be able 
to recover its ground state from sufficiently weak perturbations at sufficiently 
low temperature (the fact that we have to work with low temperatures is 
clear since we deal with the ground state). That is, we hope that there exists 
a threshold value e such that 



This condition, however, is necessary but not sufficient for error correction. It 
may happen that the m-fold degeneracy of the ground state does not survive in 
the thermodynamic limit, and we still lose some storage capacity. We therefore 
require that the ground state of the perturbed Hamiltonian remains m-fold 
degenerate for all e < eo in the thermodynamic limit. In the next section we 
elaborate further on these requirements for self-correction and show that they 
are quite easy to fulfill in a wide variety of quantum spin systems. Before we 
proceed, we must remark that these are essentially "bare-bones" desiderata; 



H A (e) :=H A + eP, 



(2) 
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e < e . 



(4) 
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more elaborate performance measures will also be investigated elsewhere [9]. 



3 Putting it together 

The main question is: which restrictions ensue on the unperturbed Hamilto- 
nian and on the perturbation PI It turns out that this question can be 
answered using the same methods that are employed for constructing low- 
temperature phase diagrams for classical spin systems with quantum pertur- 
bations [10-12]. Thus the Hamiltonian H\ can be comprised by n-spin inter- 
actions (for fixed n) that satisfy the Peierls condition [13]: the energy "cost" 
of a local perturbation u/ of a translationally invariant ground state u is of 
the order of the surface area of the region that encloses the part of the lat- 
tice on which u and uj' differ. Additionally, the unperturbed Hamiltonian is 
assumed to have a spectral gap g > (i.e., its first nonzero eigenvalue > g). 
Admissible perturbations are formed by sums of translates of an arbitrary 
Hermitian operator Pq, whose support (the set of sites on which the action of 
Po is nontrivial) is finite and encloses the origin of the lattice. Thus 



where Pe = 7(£)Po with 7(f) being the automorphism induced by the trans- 
lation of the lattice A that maps the origin to the site t and respects the 
periodic boundary conditions. Also, both the unperturbed and the perturbed 
Hamiltonians are assumed to be invariant under unitary transformations in- 
duced by a symmetry group acting transitively on the set {|cr)} fl H g A . 

Assuming these conditions are satisfied, we invoke the following theorem [10]: 
there exists a constant e such that, for all e < e , the perturbed spin system 
has m translationally invariant ground states in the thermodynamic limit. 
Furthermore, if the m translationally invariant ground states of the unper- 
turbed Hamiltonian are invariant under some additional symmetries, these 
invariance properties carry over to the ground states of the perturbed sys- 
tem. The threshold value eo of the parameter e is determined by developing a 
low-temperature expansion [14,10] of the perturbed partition function using a 
modified Lie- Trotter product formula, 



The trace is expanded in the basis {|<i)}, thus allowing for combinatorial 
analysis on a "space-time" grid, where the space axis is labelled by lattice 
sites I and the time axis is labelled by values 0, (l/N)/3, (2/N)/3, . . .,/3. The 
perturbation theory is controlled by a suitable coarse-graining of the time 
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axis, which then allows to determine the threshold value eo that will render 
the contributions of the perturbation terms Pp sufficiently small. (See [14,10] 
for details; we only remark here that translation invariance requirements can 
be lifted, the perturbation theory still goes through [10].) Similar space-time 
analysis of the error rate has been described heuristically by Dennis et al. [4]. 

Another important issue is the following: while the infinite-volume ground 
state of the perturbed system retains the degeneracy of the original (unper- 
turbed) ground state, the degeneracy may be lost in the case of finite lattice 
size. This phenomenon is referred to as "obscured symmetry breaking" [15], 
whereby low-lying eigenstates of the finite-system Hamiltonian converge to ad- 
ditional ground states in the thermodynamic limit. In this case we will have, 
for any finite A, 

AE A (e) > 0. (7) 

It is therefore important to obtain an estimate of the convergence rate in (4); 
namely, given some 5 > 0, find N such that 

AE A (e) < 5, |A| > N . (8) 

Knowing the convergence rate allows us to appraise the resources needed to 
implement error correction with the desired accuracy 5. In this respect, an 
estimate of the form 

AE A (e)=e- c ^, e < e , (9) 

where the constant c depends on e, would be ideal — an exponential gain in 
error-correction accuracy could then be achieved with polynomial resources. 
This exponential convergence rate is, in fact, one of the most attractive fea- 
tures of Kitaev's construction in Ref . [2] . 

On the other hand, the rate at which AE A (e) converges to zero is determined 
by the unperturbed Hamiltonian H A , the perturbation P, and the perturbation 
strength e. It is therefore important to know what we can expect in a generic 
setting. Obviously, the exponential behavior, as in Eq. (9), is optimal, but it 
may as well turn out that the particular implementation (e.g., with a different 
Hamiltonian) does not allow for it. We can, however, hope for a slower (but 
still quite decent) convergence rate. According to a theorem of Horsch and von 
der Linden [16], certain quantum spin systems possess low-lying eigenstates 
of the finite-lattice Hamiltonian with 

AE A (e) = c/ |A| . (10) 
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The conditions for this to hold are the following. There has to exist an order 
observable Oa of the form 

A = E°* (ii) 

where each Oe is a Hermitian operator such that [Oe, Oe>] = 0. Furthermore, 
for any interaction term $mcA in the perturbed Hamiltonian H^{e) (these 
also include the perturbation terms), we will have [$m, Oe] = unless £ G M. 
The operators $m and Oe are required to be uniformly bounded (in M and 
£ respectively), and the cardinality of the support set M must not exceed 
some fixed constant C (the latter condition has also to be fulfilled for the 
perturbation theory described above to converge). Finally, if \ip) is an eigen- 
state of H\(e), then we must have (<^|Oa|v?) = 0, but (ip\0\\ip) > (\A\ 2 (here 
the constant ( depends on Oe). The latter conditions are taken as manifesta- 
tions of "obscured symmetry breaking" [15]. Examples of systems for which 
the Horsch-von der Linden theorem holds include [15] the Ising model in the 
transverse magnetic field or the Heisenberg antiferromagnet with a Neel order. 



4 Conclusion 

Where does it all take us? It appears, from the discussion in the previous 
section, that any quantum spin system, whose Hamiltonian is formed by mu- 
tually commuting n-body interactions that satisfy the Peierls condition, can 
recover from sufficiently weak quantum (i.e., off-diagonal) perturbations at low 
temperatures. The admissible perturbations can be either finite-range [10], or 
exponentially decaying [11]. Under these (quite general) conditions, it follows 
from rigorous perturbation theory for quantum spin systems that there exists 
a critical perturbation strength eo, such that for all e < eo the degeneracy and 
symmetry properties of the ground state of the original (unperturbed) system 
survive in the thermodynamic limit. Furthermore, even if ground-state degen- 
eracy is removed by perturbation of the finite-size system, the effect of the 
error (perturbation) is effectively "washed out" in the thermodynamic limit, 
as the low-lying excited states of the perturbed system converge to additional 
ground states. 

However, the systems we have considered in this Letter are assumed to have 
classical Hamiltonians and discrete symmetries. What about truly quantum 
Hamiltonians and continuous symmetry (e.g., the quantum Heisenberg model)? 
The situation here is not so easy. For instance, it is apparent from our discus- 
sion that, in order to be self-correcting, the perturbed system must exhibit an 
order-disorder transition as the parameter e is tuned: in the "ordered phase," 
error correction is possible; in the "disordered phase," occurrence of errors re- 
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suits in irrevocable loss of information. (This has already been noted by Dennis 
et al. [4].) Since we require the ground state of the perturbed Hamiltonian to 
exhibit the same degeneracy as the corresponding state of the unperturbed 
Hamiltonian, it makes sense to talk about spontaneously broken symmetry 
in the ordered phase (i.e., when e < eo). However, according to the so-called 
Goldstone theorem [17], symmetry cannot be broken in a system with contin- 
uous symmetry and a gap. It would certainly be worthwile to explore physical 
error correction in systems with continuous symmetries as well, but the models 
in which it can work will not be as easy to find. 

Acknowledgments 

The author has benefitted from a few discussions with V. Moroz. This work 
was supported by the U.S. Army Research Office through MURI grant DAAD19- 
00-1-0177. 

References 

[I] M.A. Nielsen and I.L. Chuang, Quantum Computation and Quantum 
Information (Cambridge Univ. Press, 2000), ch. 10. 

[2] A. Kitaev, "Fault-tolerant quantum computation by anyons" (preprint 
|quant-ph/970702lD . 

[3] S. Bravyi and A. Kitaev, "Quantum codes on a lattice with boundary" (preprint 
|quant-ph/9811052| ). 

[4] E. Dennis, A. Kitaev, A. Landahl, and J. Preskill, "Topological quantum 
memory" (preprint |quant-ph/0110143 ). 

[5] A. Kitaev, Russian Math. Surveys 52, 1191 (1997). 

[6] A.L. Calderbank and P.W. Shor, Phys. Rev. A 54, 1098 (1996). 

[7] E. Knill and R. Laflamme, Phys. Rev. A 55, 900 (1997). 

[8] D. Beckman, D. Gottesman, A. Kitaev, and J. Preskill, " Measurability of 
Wilson loop operators" (preprint |hep-th/01 10205 ). 

[9] M. Raginsky, in preparation. 

[10] T. Kennedy and H. Tasaki, Commun. Math. Phys. 147, 431 (1992). 

[II] C. Borgs, R. Kotecky, and D. Ueltschi, Commun. Math. Phys. 181, 409 (1996). 
[12] N. Datta. R. Fernandez, and J. Frohlich, J. Stat. Phys. 84, 455 (1996). 



7 



[13] J. Slawny, in Phase Transitions and Critical Phenomena, vol. 11, ed. by C. 
Domb and J.L. Lebowitz (Academic Press, 1987), p. 143. 

[14] J. Ginibre, Commun. Math. Phys. 14, 205 (1969). 

[15] T. Koma and H. Tasaki, J. Stat. Phys. 76, 745 (1994). 

[16] P. Horsch and W. von der Linden, Z. Phys. B 72, 181 (1988). 

[17] L. Landau, J. Fernando-Perez, and W. F. Wreszinski, J. Stat. Phys. 26, 755 
(1981). 



8 



